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Abstract 

In the case of a two-leg Hubbard ladder we present a procedure which allows the exact deduction 
of the ground state for the four particle problem in arbitrary large lattice system, in a tractable 
manner, which involves only a reduced Hilbert space region containing the ground state. In the 
presented case, the method leads to nine analytic, linear, and coupled equations providing the 
ground state. The procedure which is applicable to few particle problems and other systems as 
well is based on an r-space representation of the wave functions and construction of symmetry 
adapted orthogonal basis wave vectors describing the Hilbert space region containing the ground 
state. Once the ground state is deduced, a complete quantum mechanical characterization of the 
studied state can be given. Since the analytic structure of the ground state becomes visible during 
the use of the method, its importance is not reduced only to the understanding of theoretical aspects 
connected to exact descriptions or potential numerical approximation scheme developments, but is 
relevant as well for a large number of potential technological application possibilities placed between 
nano-devices and quantum calculations, where the few particle behavior and deep understanding 
are important key aspects to know. 

PACS numbers: PACS No. 71.10.Fd, 71.27.+a, 73.21.-b 
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I. INTRODUCTION 



In condensed matter context, in experiments or related theoretical interpretations we 
often encounter small number of particles confined in a system or device, for example, 
in the case of quantum dots^, quantum well structures^, mesoscopic systems^, experimen- 
tal entanglement^, micro-crystals^, could gases trapped in optical lattices^, optical bound 
states^, segregationSi, interfacial stress and fracture^", self organized structures'^, sinteringi^, 
or compounds studied in the low concentration limilii^. Such problems, presenting both 
theoreticaii^*i^*i^^ and technologica l^i'°i'^i'^i'^i^° interest have continuously attracted in- 
creasing attention. Starting from even one electron problems solved exactly'^, several cases 
of interest for two2iiS2i22i24i2^^ thre o'^^i'^^1^'^1^^1^^ , fou r^°i^'i^^ , or few22i2^ particles have been 
studied, concentrating on the model behavior in the low concentration limit, or motivated 
by experimentally measured characteristics. In this hierarchy of the increasing number of 
carriers in the study of a given problem, the particle number four [Np = 4) represents a 
special case, since it is close to the particle number limit around which one can hope that 
deep rigorous descriptions can be made^S even in the non-integrable cases, the problem is 
also treatable from the numerical side as well^l, statistics and T 7^ characterization can 
be given^, and the problem retain even many-body aspects of the system's behavior— i^Siii. 

The simulations on the Np = A particle problem started more than a decade a.g(M^, but 
up today, only few valuable results are known in this subject in the condensed matter context, 
as follows. The energy dependence of the maximal Lyapunov exponent has been studied 
for ID Lenard- Jones system"^- , the spinless fermion case has been analyzed as a simplified 
model for correlated electron o^°i^' , the conjecture of the Andreev-Lifshitz supersolid has 
been studied^, entangled states have been described in the high frequency region^, doped 
quantum well structures have been investigated^', special cases where only two pairs of 
particles interact on a lattice were considered^, localization lengths have been estimated 
in ID disordered systems'^, and the behavior in the presence of Coulomb forces has been 
analyzed^. As can be seen, the knowledge accumulated in this direction is relatively poor. 
Approximated procedures have been applied in different conditions for different systems of 
interest, but the level of exact characteristics has not been reached yet. 

The need to study at exact level system holding Np = A particles is enhanced by several 
motivations. First of all, A^p = 4 it is placed in the low density limit, and as known, in this 
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limit, especially in low dimensions, no class of diagrams can be neglected in describing the 
system^^. Given by this difficulty, one often finds that traditional approximation schemes 
which work at higher densities, here fail^ or provide unphysical results^. Secondly, we 
are placed in the concentration limit where the formation of Fermi liquid properties can be 
studied^, and since this parameter region is usually numerically accessible, research with 
analytical focus, starting from numerical results, also can be done. Thirdly, as several times 
has been accentuated^!^, key aspects of the unapproximated descriptions are often hidden 
in the few particle cases. The four particle case seems to be tractable also from this point of 
view. Finally, in A^p = 4 case we face a situation which experimentally is produced, having 
potential application possibilities in several for example in the study of entangled 

states^^, non-local character of quantum theory^, high precision spectroscopy^, quantum 
communication, quantum cryptography, and quantum computation^, fields where deep and 
high quality results are clearly demanded^. 

Starting from the motivations presented above, we show in this paper that for the Np = 4 
case, exact, analytical, and explicit results holding essential information about the system 
behaviour can be indeed provided, even for arbitrary large systems. To show this, we present 
below the exact ground state for four interacting electrons placed in an arbitrary large two 
leg Hubbard ladder described by periodic boundary conditions. This is given in conditions 
for which, even the known three (quantum mechanical) particle exact results are very rare 
for systems taken outside of one dimension (see Ref. f^] and the references therein), hence 
we hope that the presented results will generate creative advancements. 

In order to obtain such results, a direct space representation is used for the wave functions. 
Starting from local particle configurations, symmetry adapted ortho-normalized basis wave 
vectors are constructed. Based on these, in the studied case, an explicit and analytic closed 
system of 9 equations is constructed, whose secular equation provides the ground state wave 
function and energy. Deducing the ground state wave function for different microscopic 
parameters of the model, ground state expectation values are calculated for different physical 
quantities of interest, and correlation functions are deduced in order to characterize the 
ground state properties. 

The method which is described here is in principle not model or particle number depen- 
dent, and could be applied for other systems as well. In presenting our calculations, the aim 
was not to hide the obtained results behind a numerical treatment or deduced symmetry 
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properties, (which certainly also can be done), but to show clear, visible, and explicit prop- 
erties which, based on the provided essential characteristics, could enhance further creative 
thinking or applications. In order to underline the importance of these aspects we note 
for example, that in recent studies made for states containing 2-4 particles, especially in 
attempts to characterize the entanglement^, or quantum dotsSS, often the analysis must be 
made without to know the state completely^i^. We show below how such ingredients, at 
least at the level of the ground state, are possible to overcome. 

The remaining part of the paper is structured as follows. Section II. presents the Hamilto- 
nian, the deduction procedure and the ground state wave functions. Section III. exemplifies 
the physical properties of the ground state. Sect. IV. presents the summary and conclu- 
sions of the paper, while the Appendices A - B presenting mathematical details, close the 
presentation. 

II. HAMILTONIAN AND GROUND STATE WAVE FUNCTIONS. 

The strategy which we use for presentation is the following one. We have chosen a simple 
model which allow us to characterize the construction of exact ground states in the presence 
of four particles. After presenting the results we indicate how the procedure could be apphed 
for other systems as well. 

A. Presentation of the Hamiltonian 

The Hamiltonian we use for presentation has the form of a standard two-leg Hubbard 
ladder Hamiltonian 

^ = -^11 E (cI,A- + ^-c-)-^± E (4A- + ^-c-) + t^E^^,Ai, (1) 

<i,j>||,o- <i,j>±,<j i 

where cj^ creates an electron at site i with spin a, ta holding the index a =||,J- are 
nearest-neighbor hopping amplitudes along and perpendicular to ladder legs, U is the on-site 
Coulomb interaction, and < i,j >a represents nearest-neighbor sites in a direction taken 
into account in the sum over the lattice sites only once. 
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B. The construction of the basis wave vectors 

If we would like to analyze by exact diagonalization the four particle problem in the singlet 
case in a two leg Hubbard ladder containing N lattice sites, we must treat numerically a 
Hilbert space of dn — [N{N — l)/2]^ dimensions, where for example at A?" = 30 we have 
dn = 2.16 • 10^, and for A?" ^ oo one encounters (ijy — >^ oo as dj? ~ A^^. 

We show below how it is possible to deduce exactly the ground state for a such type of 
system in the case of an arbitrary large two leg Hubbard ladder based on only nine linear 
and analytic equations, and to extract essential information from the obtained results. In 
order to do this, first of all we dehmit exactly the Hilbert space region (Hg) containing the 
ground state by the construction of nine type of orthogonal basis wave vectors spanning Tig. 
This procedure is presented below. 

1. The generating configurations 

We are interested first to have an image about the possible type of states of the studied 
four particles in the system under consideration. To obtain such type of information, we 
number all lattice sites of the ladder as shown in Fig.l (periodic boundary conditions are 
considered). In the figure, A^, considered even number, denotes the number of sites within 
the system, while n = N/2 gives the number of rungs, respectively. Using now an r-space 
representation, one observes that since the ladder legs, and the spin reversed configurations 
are equivalent, the studied four particles can be placed into the system only in nine possible 
ways, as depicted in Fig. 2. The presented possibilities, denoted by capital letters A to J will 
provide nine type of basis wave vectors (denoted by the same letters), whose construction 
is presented below. We mention that the subscripts i,j,k are denoting particle positions 
within the considered states A to J presented in Fig. 2, which are such chosen, to have the 
first particle position placed into the origin (e.g. lattice site 1). In the following, the nine 
possible four-particle states presented in Fig. 2 will be called generating configurations. How 
one arrives from the generating configuration X — A,B,...,J to the base vector \X), is 
explained in the following two subsections. 
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2. The sum of configurations related to each generating configuration 

Each generating configuration provides other seven related configurations (brother con- 
figurations) of the same type. These are obtained by a) rotating the generating configuration 
by 180 degrees along the longitudinal symmetry axis of the ladder, b) rotating the gener- 
ating configuration by 180 degrees along the symmetry axis perpendicular to the ladder, c) 
rotating by 180 degrees the configuration obtained at b) along the longitudinal symmetry 
axis of the ladder, and finally, d) other four related configurations are obtained by revers- 
ing all spin orientations in the generating configuration and the configurations deduced at 
points a)-c). As an example, the eight related configurations describing the state Di j taken 
at i = 2, j = 3, are depicted in the first column of Fig. 3. 

After this step, since all lattice sites are equivalent, the different ,, related" configurations 
are translated by elementary translation N/2 times along the ladder, and all the contri- 
butions are added. We obtain in this manner a sum of configurations for each generating 
configuration. Such a sum contains 8 x N/2 components. For example, in the case of the 
1)2,3 state, this sum is presented in Fig.3. 

The procedure described above must be effectuated separately for each generating con- 
figuration. As a result, we obtain at this point nine configuration sums. Each of these sums 
will give rise to one basis wave vector as follows. 

3. The basis wave vectors 

A given configuration sum described in the previous subsection provides one basis wave 
vector if each individual configuration of the sum is written in mathematical form via four 
creation operators acting on the bare vacuum. In order to do this, we have to fix the order 
of creation operators for each type of contribution, which has been done as follows. For two 
doubly occupied sites we write the creation operators of the couples next to each other, first 
the spin up, then the spin down contribution, as c]|cj j^c] .|.c] JO), where only the restriction 
i ^ j exists. In the case of basis wave vectors containing only one doubly occupied site at i 
one uses cj |c| j^ct -i-c^ JO), were i j and i ^ k must hold. Finally, for cases without double 
occupancies, the convention c| |c] ^c^ j^cj JO) is considered, where i < j, and k < I must hold. 
Using these conventions, for example, in the case of |-Djj), taken ai i — 2, j ~ 3 and depicted 
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in Fig. 3. the result becomes 

1-02,3) = ((cI|cJ^4TC(n+3)i + 4T4i4T^{n+4)i + '^h'^ll'^h'^ln+5)l + • • •) 
+ (4|4ic|„4.3)|4i + 4T4i^|„+4)|C3i + 4T4iC(„+5)|C4i + • • •) 

+ (C(n+l)T4"+l)i^3T4«+2)i + 4n+2)T4"+2)i^4T4n+3)i + 4n+3)T^(n+3)i4TC(n+4)i + • • •) 

+(4|4i4T4n+i)i 4T4i4T4n+2)i + 4T4i4T4n+3)i + • • •) 
+(4T4i4"+i)T^2i + 4T4i4"+2)T^3i + 4T4i.4"+3)T^4i + • • •) 

+ (4"+3)T^("+3)i4"+2)T^U + ^(n+4)T^(n+4)i4n+3)T^2J. + 4"+5)T^(n+5)i4«+4)T^3i + • • •) 

+ (^(n+3)T^(n+3)i^lT^(n+2)i + C(„+4)|4n+4)i4T^(n+3)i + ^(n+5)T^(n+5)i4T^(n+4)i + ' ' •))|0) 

Similar procedure applies for all basis wave vectors. We mention that the so obtained basis 
wave functions are orthogonal. 

Here we must note that because of the fixed conventions presented above, sometimes 
an additional negative sign arises in the process of writing the mathematical expression 
corresponding to a basis wave vector component translated from the end to the beginning of 
the ladder in the presence of the periodic boundary conditions. For example, if we translate 
the vector 4,t^!v/2 t^2,j,4,j,l^) elementary translation along the ladder, according to 

the fixed conventions one obtains 4,t4,t4,j.^4,il'-') ~ ~4,t4,t4,j.^4,J.I^)' 



C. The ground state wave function 

After the calculation presented above, we are in the possession of nine type of orthogonal 
basis wave vectors \Ai), \Bi), | Jij,fc), enumerated together with their generating configu- 
ration in Fig.2. Let as denote these basis wave vectors by ), m = 1, 2, 3, 9. Now one 
observes that by applying the Hamiltonian on a given \ol™'j) basis wave vector with fixed 
m, we obtain the result inside the )} set. Consequently, nine explicitly given analytic 

linear equations form a closed system of equations, whose secular equation, by its minimum 
eigenvalue, contains the ground state at attractive U. The nine equations are exemplified in 
Appendix A and are available in their complete extent in Ref. [^^]. The ground state nature 
of the minimum energy eigenstate has been tested by exact numerical diagonalizations taken 
on the full Hilbert space for different values. 
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The fact that the analytic solution of the problem can be given in such a manner for 
arbitrary large ladder length is connected to the observation that with increasing A^, the 
type of the particle configurations describing the system (see Fig. 2), remains unchanged. The 
deduction of the ground state itself from the system of equations presented in Appendix A 
must be numerically given^. Since the possible inter-particle distances (e.g. the possible 
values of the indices in O-^^ at fixed m) depend on the A^ value, the number of 

equations which must be numerically treated depends on A^ in the frame of the same analytic 
expressions. For example, for the m = 1 case we have 1 < i < 1 + N/A, for the m = 2 case 
we have 1 < < 1 + N/A, etc. The number of obtained equations de is however significantly 
lower than dn, the Cg = dn/de ratio being at least of order 10^ at intermediate A^ ~ 0(10) 
values. Increasing A^, Cg further increases. 

D. Application possibilities in other cases 

In fact, the deduced system of equations, based on symmetry properties, delimitates 
from the full Hilbert space a d^, dimensional space region, inside of which the ground state 
is placed. The deduction of a such region is possible for other (non-disordered) models, and 
other particle numbers as well. In order to do this, we mention that if the lattice sites are 
equivalent, the elementary translation of a particle configuration can be in principle given 
with a site independent multiplicative phase factor exp{iatrans) ■ Furthermore, the rotation of 
a particle configuration along a symmetry axis can be given in principle with a multiplicative 
phase factor of the form exp{iarot), both atrans, c^rot providing their contributions in the basis 
wave vectorsSi. In the described case, we have atrans = Oirot = 0, but in other cases, the 
energy can be minimized in function of these parameters. 

In deducing the linear system of equations describing Tig in a new case characterized by a 
new H, one must start from a given basis wave vector (denoted by for example). This is 
obtained from a generating particle configuration, which is translated and rotated as specified 
above, all such obtained configurations being summed up. From technical reasons, the first 
generating particle configuration must be such chosen to contain (for 1/2 spin fermions) only 
double occupied sites placed in nearest neighbor sites. Calculating now H\vi), the result 
becomes a linear combination containing new base vectors |f2),..., |fni), holding the same 
symmetry properties, but being related to new generating configurations. Continuing the 
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procedure by calculating H\v2), \H\v-i), etc., since periodic boundary conditions are used, 
the linear system of equations closes up. It is even not important to know all distinct particle 
configuration possibilities, since these are automatically generated by the H\vi) operation. 

III. GROUND STATE PROPERTIES 

By diagonalizing the system of equations presented in Appendix A and taking the mini- 
mum energy solution, one finds the ground state wave function I^E^^). Using this, the complete 
quantum mechanical characterization of the ground state can be given. In order to exem- 
plify the results, we present in (IB1IB2|) explicit expressions containing the leading terms of 
the ground state wave function for two parameter values. Even Appendix B shows that in 
the leading terms of the ground state wave function, the particles have the tendency to be 
placed in pairs, the pairs tending to occupy the highest possible distance between them. 
This is reflected as well in the density-density correlation function depicted in Fig.4c. 

Ground state expectation values and correlation functions are exemplified in Figs. 4-5. 
calculated for = 28, e.g. ladder containing 14 rungs described by periodic boundary 
conditions taken along the ladder. The correlation functions are defined as follows. The 
density-density correlation function has the expression 

1 ^ 

Cn{r) = — Y,{{nini+r) - {ni){hi+r)) (2) 

i=l 

where -hi = + fii^cr = c\^Cia- The spin correlations are studied via 

1 ^ . . , . 

Cs^{r) = J2iiS!SUr) - {S!){SUr)), (3) 
i=i 

where = (l/2)(nj,^ — hi^i). The superconducting pairing s-wave^ correlation function is 

1 ^ 

CsUr) = -TT Y.i(4AAi+r)Ai+r)'\) " (clTC(i+r)T ) (cliC(i+r)i) ) , (4) 
i=l 

and the superconducting pairing d-wave^'^ correlations are studied via 

1 ^ . 

C,Ur) = ^T.(^\^ + r)A{^) (5) 
1=1 

where A{i) = (cj^iCj^-f — Cj2|Cjj|). The i in A{i) denotes a rung connecting the lattice sites 
ii, 12- The r values inside the figures are given in lattice constant units. 
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Fig.4a presents the ground state energy and the potential energy in t\\ units in function 
oi u = \U/t\\\ at t|| = t±. Fig. 4b shows that the spin-spin correlations are exponentially 
decreasing, the decrease rate in the exp{—r/^) being of the form 1/^ = 0.34 + 0.78^/jitj. 
The density- density correlations depicted in Fig.4c show that the particles tend to occupy 
opposite positions in the ladder closed by periodic boundary conditions. 

In Fig. 5 the behavior of the superconducting correlation functions is presented. In these 
plots u = U/t\\ holds. The correlations in Fig. 5 are decreasing with r, and for s-wave 
case slightly increase by increasing the attractive U value. Fig. 5c further shows that the 
decrease of the inter-leg hopping amplitude at fixed on-site interaction is detrimental to 
d-wave pairing correlations. Similar behavior has been found also by others^. 



IV. SUMMARY AND CONCLUSIONS 



We describe a procedure which allows the exact deduction of ground state wave functions 
for few particles in lattice models. The main result of our paper is that indeed, a such type 
of analytic description can be done. In the case of an arbitrary large two leg Hubbard ladder 
taken with periodic boundary conditions and containing four electrons, presented in details, 
the method leads for the singlet state to nine analytic linear and coupled closed system of 
equations, whose secular equation, through its minimum eigenvalue solution, provides the 
ground state wave function and ground state energy. The procedure is based on a r-space 
representation of the wave functions and properly constructed symmetry adapted orthogonal 
basis wave vectors. These are obtained from generating particle configurations translated 
and rotated in the lattice and finally added. The linear system of equations is obtained by 
applying the Hamiltonian on the deduced basis wave vectors. The procedure can be applied 
for other systems as well. 

The fact that the analytic structure of the ground state becomes visible by the use of 
the method underlines that the presented procedure contributes not only to the understand- 
ing of theoretical aspects related to exact descriptions, or development possibilities of new 
numerical approximation schemes, but has implications on a broad spectrum of subfields 
related to technological developments placed in between nano-devices and quantum compu- 
tation, where the exact knowledge of the behavior of a small number of quantum mechanical 
particles plays a main role. 
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APPENDIX A: THE LINEAR SYSTEM OF EQUATIONS CONTAINING THE 
GROUND STATE. 



This Appendix presents the nine analytic equations describing the action of the Hamil- 
tonian on the basis wave vectors. 

The first two equations are devoted to the \Ai), \Bi) species containing only (two) doubly 
occupied sites. 

H\Ai) = 2u\Ai) — t±\Di^i) — Ii>2\Ci-i^i) — /i<||Ci,i+i) , 
H\Bi) = 2u\Bi) - t±Ii^i\Di^i) - Ii^i\Ei_i^i) - /j<n|£^j_j+i) , 

where //^ = 1 if the statement K is true, and Ik = otherwise. 

The following three equations describe the action of H on the basis wave vectors contain- 
ing only one doubly occupied site (|Cij), \Dij), \Eij)) as follows 

H\Cij) = 

- 4(5,-,+i|A) - 4(5,-,+i(l + - (1 - k2)\Ci-i,j) 

-(1 - - ^i,2Sj,ii±i+i)\Ci,j-l) 

~(1 — <^«,2<^,- H+l + Sj,n-i+l — Sj,n-i+2)\Ci,j+l) 
~(1 ~ + <^j,n-i+l - (^j,n-j+2)|C'j+lj) 

+5j,2(l + (^j,3)(l - (^j,2±i - (^j,2±l+i)|C'2,n-j+3) 

-tJA,,)-*x 



-(1 - ^1,2) (1 + 5j,i+l)|-^i-l,i,n-j+i+l) 



+ 
-( 



^ (1 - 5,,2)(1 + W)(l - M,>^)l^^,2j) 1 

— Si :ilj^ n+i+l 9,.n,-j+i+l) J 

(1 + + Sj^n-i+1 - ^j,n-i+2)\Fi,i+l,n-j+i+l) 



11 



+ (1 + + 5j^n-i+l - (^j,n-i+2) X 

(1 - (^i,2-^j>Il±i)|-^i+l,2j+l) 
~5i,2-^j>2±i |-fi+l,2,n-j+i+l) 



X < 



— 2 

-^j>Il+i+i |G^i,n-i+i+l,i) 



-^±(1 — <^j,n-i+2) ■ ' 



^i<f +1 ■ 



^j<2i-l 



/j>2i-l|Gj,i,l) 



^j<2i 



^i>2i-l 



-l|(jn- 



-j+2,n-i+2,l/ 



■i+2,n-i+i+l,n-i+2 



-t±{l - Sj^i - 5j^i + 5j,n-i+2) ■ < 



+ < 



Ij<n-i+2lj<i{\Cj^i) + \Ej^ij) 

In-i+2>j>i 

Ij>n-i+2{\Cn-j+2,n-i+2) + |-E'n-i+2,n-i+2)) 

^i>l[l - '5i,2(5i,2 + ^i,Il±i+l) + ^i,f+l5i,2]|A,i-l) 
- 5i,f+l)|A,n) 

[1 - Si,2{Sj,l + Sj n±i) - 5j,n + Sin+l{SjR - SjR+l)]\DiJ+l) 

(<^i,t-^j>f +1 + ^i,t+l(l - - A-i+l,n-j+2) 

[1 - (^i,2 - kfi{^j,2 + (^j,3 + /j>S+i)]|Gi,2,i) 

— ^i,3-(?>f +2|Gj^2,n-i+i+l) 



— < 



-(1 - C>i,2) { 

\ -(7>«l^i-l,i,n-j+i+l) 



+ < 



(l-<^.f+l) 



^i+l.i-j+iy 



^ (^i,f+l-^l<j<f+l(l - (^n,4(^j,2)|Gi,2,n-j+2) ^ 
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-5,: 



i,2 



-t± ■ { 



[1 - 6j^n - +1 - Si,2{Sj,l + 5j,2 + Ij>^+l)]\G,+i^2,, 
'5i,n.|G'j+l,2,l) 

2 + l/l<j<£ + l|G's S+l S+j) 

^^■,21^4+1, 2, l) 

-^f +l<j<n|G'j+i_2,n+j-j+l) 
h<j<i\Hjj>n-i+j+l) 

i,n-j+i+l) 
Ij>n-i+2\Hn-j+2,n-j+2,n+i-j+l) 



(1 ~ ^j,i){h<j<n-i+2 + 25j^n_.j+2)|^j,lj) 



^j>n-i+2 



\Jn-j+2,l,n-i+2) 



H\E^,) = 

u\E,^,) - + + (1 + 5. - t^{l - 

[1 - 6i^i + 5i,2('5j,f+i - Ij>^+i)]\Ei_ij) 



+ 



/j<n+l|L'j-i) 



-^f +l<i<n-i+2|-D„_j+2,n-i+2) 



~ ' Si,2lj>^+l\Ei-i^n~j+2) 

^ + 5j,2 - 5j^^+i)\E2,n-j+2) 

^ l^i.j+OI-SiJ-l) ^ (1 + — Sj,n-i+2 + " <^i,l<^j,f +l) 

— / n+i + l 
- 2 



+ (1 + Sj,i+i] 



I -^i<2i-l|Gjj_,;+ij) 
^Ij>2i-l\Gj^i^l) 

> 

—Ij<2i-l\Gn~j+2,n-i+2,l) 
-'i>f +1 ■ \ } 
[ -^j>2i-l |G'„_j+2,n-j+i+l,n-j+2) J 

-^i>2| A-l,n-j+i+l,i) 
'^i,2|-f^l,2,n-i+i+l) 

25j,2|-f^2,2,l) 



In- 



-l>j>l|-f^2,2 



^ + [l + 5,-,+i(l + 25,i)]l^^.2,,) 



+ l^i.i+l — +l)|-f^l,2,n-i+j+l) 

-5j-„_,+2)|ifm,2,i+i). 

The last four equations devoted to the base vectors l-F^j./c), \Gij^k), \Hij^k), \Ji,j,k) (not 
containing doubly occupied sites) can be find in Ref.p^i]. 
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APPENDIX B: EXEMPLIFICATION FOR GROUND STATE WAVE FUNC- 
TIONS 



We present below the leading terms of explicit ground state wave functions deduced 
for = 28, at \U/t\\\ — 3. The ground state \'^g) is normahzed to unity, and contains 
ortho-normalizcd basis wave vectors. 

For = 0.8 one obtains for the ground state wave function 

1^.) = 

0.181883|£;7,8) + 0.181878|C7_8) + Q.n^7m\DTj) + 0. 169247] Cgj) 
+0.169246|£;6,7) + 0.157289|D6,6) + 0.145021|C5,6) + 0. 145004 1£;5,6) 
+0.138346|L'8,7) + 0.1383461^7,8) + 0.128723|L'6,7) + 0.1287211^7,6) 
+0.12823|L'5,5) + 0.1113461^4,5) + 0.111315|£;4,5) + 0.110256|L'5,6) 
+0.1102391^6,5) + 0.10177|£;6,8) + 0.101761|C6,8) + 0.097877] ^7,8,1) 
-0.0978768]G8,2,8) + 0.0931102]L)7,9) + 0.0913856]L'4,4) - 0.0910721]G7,2,7) 
+0.0910714]G'6,7,i) + 0.0909671 ]C5,7) + 0.0909645 ]£;5,7) + 0.0898481 ]i:'8,6) 
+0.0898458] ^6,8) + 0.0845566]D4,5) + 0.0844776 ]D5,4) + 0.0802931]D5,7) 
+..., (Bl) 

while for t±/t\\ — 0.1 one has 

0.29866]E7,8) + 0.29365]E6,7) + 0.284039]^5,6) + 0.270726]E4,5) 
+0.255311]£;3,4) + 0.240496]£;2,3) + 0.230556]Ei,2) + 0.16658]C7,8) 
+0.156952]C6,7) + 0. 149053 ]£;6,8) + 0.145317]£;5,7) + 0.13942]£;4,6) 
+0. 137923 ]C5,6) + 0.131927]E3,5) + 0. 123835 ]E2,4) + 0.116908]Ei,3) 
+0.110056]C4,5) + 0.081715]C6,8) + 0.0755983]E6,9) + 0.0751106]i/6,i3,7) 
-0.0750843]//7,2,9) + 0.0746358]C3,4) + 0.0744796] (^5,7) + 0.0743108] £;5,8) 
+0.0740908] 57) - 0.073841 ]i/5,i3,6) - 0.0738214]i/6.2,8) - 0.0722444]B6) 
+0.0718117]E4,7) - 0.0713933]i/5,2,7) - 0.0713921]if4,i3,5) + 0.0693442]55) 
+... (B2) 
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FIG. 1: The numbering of the lattice sites for the two leg ladder taken with periodic boundary 
conditions. N denoting the number of lattice sites is considered even. The t±, (t||), denotes the 
inter- leg, (intra- leg) hopping matrix element. 
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FIG. 2: The different possible types of base vectors. We note that for the cases C,E i ^ j, while for 
F, J j < k is considered, respectively. In the cases F, G, H, J, the double occupancy is forbidden. 
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FIG. 3: The structure of the |-D2,3) base vector. 
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Figure caption for Fig.4: 

The properties of the ground state for t± — t\\. (a) The dependence of the energy (in 
t\i units) on M = f//t||. The continuous hne is the total energy, while the dots indicate the 
potential energy, (b) The logarithm of the same-leg S^-S^ correlation function for -u = 
(dots, dot-dashed line), -u = — 10 (squares, long dashed line), u = —30 (diamonds, short 
dashed line), u — —100 (stars, continuous line), (c) The same-leg density-density correlation 
function ior u — (dots, dot-dashed line), u — —10 (squares, long dashed hne), u — —30 
(diamonds, short dashed line), u — —100 (stars, continuous line). 
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Figure caption for Fig. 5: 

Superconducting ground state correlation functions, (a) The same-leg superconducting 
s-wave correlation function for t± — i|| and u — (dots, dot-dashed line), u — —10 (squares, 
long dashed line), u = —30 (diamonds, short dashed line), u = —100 (stars, continuous 
line), (b) The superconducting d-wave correlation function for t± = t\\ and u = (dots, 
dot-dashed line), u — —10 (squares, long dashed line), u — —30 (diamonds, short dashed 
hue), u — —100 (stars, continuous hue). We mention that the curves corresponding to the 
last two u values are almost superposed, (c) Superconducting d-wave correlation function 
for u = —10 and t = t_L/t\\ taken ast — 1 (squares, dot-dashed line), t — 0.5 (triangles, short 
dashed line), t = 0.3 (X-s, long dashed line), t = 0.01 (circles, continuous line). 
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